2101.03179v1 [astro-ph.GA] 8 Jan 2021 


arXiv 


DRAFT VERSION JANUARY 12, 2021 
Typeset using IATEX twocolumn style in AASTeX63 


A Luminous Quasar at Redshift 7.642 


FEIGE WANG,!* JiNvi YaNaG, ^! Xraonur FAN,! Josepu F. HENNAWL?? AARON J. BARTH,* EDUARDO BANADOS,? 
FuYAN BIAN,’ KONSTANTINA BoursiA,Ó THOMAS CONNOR,’ FREDERICK B. DAVIES,’ ROBERTO DECARLI,® 
ANNA-CHRISTINA EiLERS,: * EMANUELE PAOLO FaRINA, RICHARD GREEN,! LINHUA JIANG,” JIANG-Tao Li," 
CHIARA MAZZUCCHELLI,” RICCARDO NANNL? JAN-TORGE SCHINDLER,? BRAM VENEMANS,? FABIAN WALTER,’ 
XvuE-BiNa. Wu,!?^!^ AND MiNGHAO YUE! 


1 Steward Observatory, University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721, USA 
? Department of Physics, University of California, Santa Barbara, CA 93106-9530, USA 
3 Maz Planck Institut für Astronomie, Königstuhl 17, D-69117, Heidelberg, Germany 
^ Department of Physics and Astronomy, University of California, Irvine, CA 92697, USA 
5 European Southern Observatory, Alonso de Córdova 3107, Casilla 19001, Vitacura, Santiago 19, Chile 
6 Las Campanas Observatory, Carnegie Observatories, Colina El Pino, Casilla 601, La Serena, Chile 
T Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, Pasadena, CA 91109, USA 
8 Lawrence Berkeley National Laboratory, 1 Cyclotron Rd, Berkeley, CA 94720, USA 
9 INAF — Osservatorio di Astrofisica e Scienza dello Spazio di Bologna, via Gobetti 93/3, I-40129 Bologna, Italy 


10 MIT-Kawli Institute for Astrophysics and Space Research, 77 Massachusetts Avenue, Building 37, Room 664L, Cambridge, 
Massachusetts 02139, USA 


11 Maa Planck Institut für Astrophysik, Karl-Schwarzschild-Strafie 1, D-85748, Garching bei München, Germany 
1? Kavli Institute for Astronomy and Astrophysics, Peking University, Beijing 100871, China 
13 Department of Astronomy, University of Michigan, 311 West Hall, 1085 S. University Ave, Ann Arbor, MI, 48109-1107, USA 
14 Department of Astronomy, School of Physics, Peking University, Beijing 100871, China 


ABSTRACT 


Distant quasars are unique tracers to study the formation of the earliest supermassive black holes 
(SMBHs) and the history of cosmic reionization. Despite extensive efforts, only two quasars have been 
found at z > 7.5, due to a combination of their low spatial density and the high contamination rate in 
quasar selection. We report the discovery of a luminous quasar at z — 7.642, J0313—1806, the most 
distant quasar yet known. This quasar has a bolometric luminosity of 3.6 x 10? Lo. Deep spectroscopic 
observations reveal a SMBH with a mass of (1.6+0.4) x 10? Mc in this quasar. The existence of such a 
massive SMBH just ~670 million years after the Big Bang challenges significantly theoretical models of 
SMBH growth. In addition, the quasar spectrum exhibits strong broad absorption line (BAL) features 
in C IV and Si IV, with a maximum velocity close to 20% of the speed of light. The relativistic BAL 
features, combined with a strongly blueshifted C IV emission line, indicate that there is a strong active 
galactic nucleus (AGN) driven outflow in this system. ALMA observations detect the dust continuum 
and [C II] emission from the quasar host galaxy, yielding an accurate redshift of 7.6423 + 0.0013 and 
suggesting that the quasar is hosted by an intensely star-forming galaxy, with a star formation rate 
of ~ 200 Mc yr^! and a dust mass of ~ 7 x 10" Mọ. Followup observations of this reionization-era 
BAL quasar will provide a powerful probe of the effects of AGN feedback on the growth of the earliest 
massive galaxies. 


1. INTRODUCTION 


Luminous high-redshift quasars are key probes of the 
history of cosmic reionization. Deep spectroscopy of 
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z > 6 quasars indicates that the intergalactic medium 
(IGM) is significantly neutral at z = 7 (e.g. Banados et 
al. 2018; Davies et al. 2018; Wang et al. 2020a; Yang 
et al. 2020a) but highly ionized at z S 6 (e.g. Yang 
et al. 2020b), providing crucial information to map the 
cosmic reionization history in a way complementary to 
the cosmic microwave background (CMB) integral con- 
straint on reionization from measurements of the elec- 
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tron scattering optical depth (e.g. Planck Collaboration 
et al. 2020). 

In addition, the earliest supermassive black holes 
(SMBHs), the engines of the most distant quasars, are 
crucial for understanding the formation mechanisms of 
the first generation black hole seeds (see Inayoshi et al. 
2020, for a recent review). A billion solar mass SMBH 
at z ~ 7, having grown at the Eddington limit since for- 
mation, requires a seed black hole of mass ~ 1000 Mo at 
the time the first luminous object formed in the Universe 
(i.e., z ~ 30, Tegmark et al. 1997). This growth corre- 
sponds to a factor of ~ 10 increase in mass within a 
mere ~ 650 Myr. The recent discovery of a 1.5 x 10? M 
SMBH in a luminous quasar at z — 1.52 poses the most 
stringent constraints yet on the masses of the seed black 
holes (Yang et al. 20202). 

After a decade of industrious searches, a sample of 
more than 50 quasars now exists at 6.5 « z « 7 (e.g. 
Mazzucchelli et al. 2017; Wang et al. 2019; Yang et al. 
2019; Reed et al. 2019; Matsuoka et al. 2019a), enabling 
the first statistical studies of the early accreting SMBHs 
deep into the epoch of reionization (Wang et al. 2019). 
However, the sample of quasars at z > 7 is still very 
limited, since the Lya emission is redshifted to near- 
infrared wavelengths, making both imaging and spec- 
troscopic observations more challenging. In addition, 
the number density of z > 7 quasars is low (Wang et al. 
2019) while the contaminants, Galactic cool dwarfs and 
early type galaxies, are far more numerous. To date, 
there are only seven quasars known at z > 7 (Mortlock 
et al. 2011; Banados et al. 2018; Wang et al. 2018; Yang 
et al. 2019, 2020a; Matsuoka et al. 2019a,b) with two 
of them at z = 7.5 (Barados et al. 2018; Yang et al. 
20202). 

In this Letter, we report the discovery of a new 
quasar J031343.84—180636.4 (hereinafter J0313—1806) 
at z — 7.642, the most distant quasar known to date. 
The existence of a 1.6 x 10? Mc SMBH and strong out- 
flows in this quasar provide new insights into the forma- 
tion and the growth of the earliest SMBH. Throughout 
this Letter, we use a ACDM cosmological model with 
Ho = 70.0 km s^! Mpc^!, Qm = 0.3, and Q4 = 0.7 and 
the AB photometric system. 


2. OBSERVATIONS AND DATA REDUCTION 


J0313-1806 was selected as a z > 7.2 quasar candidate 
from our ongoing reionization-era quasar survey (Wang 
et al. 2017, 2018, 2019; Yang et al. 2019, 2020a). This 
survey relies on an imaging dataset that combines opti- 
cal imaging from Pan-STARRSI (PS1, Chambers et al. 
2016) and the DESI Legacy Imaging Surveys (DELS, 
Dey et al. 2019), as well as infrared imaging from the 


UKIRT Hemisphere Survey (UHS, Dye et al. 2018), 
the VISTA Hemisphere Survey (VHS, McMahon et al. 
2013), and the Wide-field Infrared Survey Explorer sur- 
vey (WISE, Wright et al. 2010). J0313-1806 falls into 
a sky area covered by PS1, DELS, VHS, and WISE. 
It drops out in all of the optical bands but is well de- 
tected in the infrared bands (see Figure. 1) with colors 
of z—J > 3.7, J — W1 = 0.91, and W1 — W2 = —0.21. 
'The sharp dropout in the z-band and the blue slope of 
the continuum make it a promising quasar candidate at 
z > 7.2. The detailed photometry for J0313-18006 is 
listed in Table 1 and the cutout images are shown in 
Figure 1. 


2.1. Near-infrared Spectroscopy 


'The initial spectroscopic observation for J0313-1806 
was obtained on 2019 November 04 (UT) with Mag- 
ellan/FIRE (Simcoe et al. 2010) using the high- 
throughput longslit mode. A 15-minute exposure shows 
a clear Lyman break at 1.048 um, indicating that it is a 
source at z = 7.6. We then followed this object up with 
Magellan/FIRE in Echelle mode and with the JH grism 
on Gemini/Flamingos-2 (Eikenberry et al. 2004) and 
confirmed it as a high redshift quasar. Since the Mg II 
emission line, the most reliable line for measuring black 
hole mass at high redshifts, is redshifted to the edge of 
the ground-based infrared observation window, exten- 
sive spectroscopic observations were obtained with Mag- 
ellan/FIRE (Echelle mode), Gemini/Flamingos-2 (with 
HK Grism), Keck/NIRES! and Gemini/GNIRS (Elias 
et al. 2006) to improve the S/N at this wavelength. De- 
tailed information for all the observations is listed in 
Table 1. 

All spectra were reduced with the spectroscopic data 
reduction pipeline PyPeIt? (Prochaska et al. 2020a,b). 
The pipeline includes flat fielding, wavelength calibra- 
tion, sky subtraction, optimal extraction, flux calibra- 
tion and telluric correction. More detailed descriptions 
of each processing step can be found in Prochaska et al. 
(2020b) and Wang et al. (2020a). We stacked all spec- 
tra obtained from FIRE (Echelle mode), Flamingos-2 
(with the HK Grism), NIRES and GNIRS after binning 
them to a common wavelength grid with a pixel size of 
~ 86 km s^! (similar to the GNIRS native pixel scale). 
We then scaled the stacked spectrum by carrying out 
synthetic photometry on the spectrum using the VISTA 
J-band filter response curve to match the J-band pho- 
tometry for absolute flux calibration. The final spec- 
trum after correcting for Galactic extinction based on 


1 https://www2.keck.hawaii.edu/inst /nires/ 
2 https: //github.com/pypeit /Pypelt 
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Figure 1. Upper panel: image cutouts (20” x 20", north is up and east is to the left) for J0313-1806 in PS1 z, PSL y, 
DELS z, VISTA J, VISTA Ks, WISE W1 and WISE W2 bands. The photometry is given in Table 1. Lower panel: the final 
stacked spectrum of J0313-1806. In the figure, we re-binned the spectrum by two spectral pixels (~ 173 km s^!) for illustration 
purposes. The black and gray lines represent the Galactic extinction-corrected spectrum and the error vector, respectively. The 
blue line denotes the quasar composite spectrum constructed with SDSS z ~ 2 quasars having similar C IV blueshifts and line 
strengths. The purple dashed line denotes the power-law continuum. The orange points are flux densities determined from 
photometry in the J and Ks-bands. The inset panel shows the Mg II line fitting with the purple dot-dashed line denoting 
the power-law continuum, the green dashed line denoting the pseudo-continuum model (the sum of power law continuum, Fe II 
emission, and Balmer continuum), the orange line representing the Gaussian fitting of the Mg II line and the red line representing 
the total fit of pseudo-continuum and Mg II line. The thin grey lines in the insert panel represent the spectral fitting of 100 
mock spectra as described in 83. 


the dust map of Schlegel et al. (1998) and extinction 
law of Cardelli et al. (1989) is shown in Figure 1. 


2.2. ALMA Observations 


In order to determine the systemic redshift and inves- 
tigate the host galaxy properties, we observed J0313- 
1806 with the Atacama Large Millimeter/submillimeter 
Array (ALMA) in the C43-4 configuration (program ID: 
2019.A.00017.S, PI: F. Wang). We tuned two spectral 
windows centered at the expected frequency of the [C II] 
line and the other two spectral windows centered at 
about 15 GHz away from the expected [C II] line. The 
observations were taken in 2020 March 02 with 29 min- 
utes of on-source integration time. 

The ALMA data were reduced using the CASA 5.6.1 
pipeline (McMullin et al. 2007) following the standard 
calibration procedures. The continuum map was im- 
aged by averaging over all the line-free channels using 
Briggs cleaning via the CASA task tclean with robust- 


ness parameter r — 2.0, corresponding to natural visibil- 
ity weights, to maximize the signal-to-noise ratio (S/N) 
of our observations. The beam size for the continuum 
image is 0770 x 0756. We subtracted the continuum us- 
ing the uvcontsub before imaging the [C II] line. The 
(C I] map was collapsed over c1.4ejgye to maximize 
the S/N of the intensity map. Since both the contin- 
uum and [C II] emissions are marginally resolved, we 
extracted the 1D spectrum with an 1"5 diameter aper- 
ture centered on the emission. The fully calibrated con- 
tinuum map, [C II] map and [C II|2-continuum spectrum 
are shown in Figure 2. From the spectral fitting of the 
[C II]--continuum spectrum, we derive a [C II] redshift 
of zic = 7.6423 + 0.0013. 


3. A 1.6 BILLION SOLAR MASS BLACK HOLE 


The most reliable tool for measuring the mass of 
SMBHs at high redshift is the Mg II virial estimator 
(Vestergaard & Osmer 2009). We first de-redshift the 
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Figure 2. ALMA observations of the dust continuum and [C II] line. The left panel shows the dust continuum map with 
contour levels of [3, 5, 7, 9, 11, 13] xo, where ø = 0.02 mJy. The middle panel shows the integrated [C II] emission with contour 
levels of [3, 4] xo, where ø = 0.04 Jy km s^!. The sizes of the continuum cutout and the [C II] cutout are 4" x 4". The [C II] 
intensity map was collapsed over +1.40jine. The white crosses in both the left and middle panels highlight the quasar position 
derived from VHS infrared images. The right panel shows the [C II] spectrum (black) and noise vector (blue) extracted from the 
data cube with an aperture diameter of 1"5 centered at the peak position of the continuum map. The spectrum is re-binned to 
23.4 MHz channels (~ 30 km s~'). The spectral fitting (red line) gives FWHM = 312 +94 km s^! and zjcr = 7.6423 + 0.0013. 


spectrum to the rest-frame using the [C II] redshift. 
Then we fit à pseudo-continuum model that contains 
a power-law continuum, Fe II emission (Vestergaard & 
Wilkes 2001; Tsuzuki et al. 2006) and a Balmer contin- 
uum (De Rosa et al. 2014) to spectral regions free of 
strong, broad emission/absorption lines (except for the 
Fe II). This procedure allows us to measure a rest-frame 
UV slope of a, = —0.91 + 0.02 and a quasar bolometric 
luminosity of (1.4 + 0.1) x 104" erg s~! after applying a 
bolometric correction factor of C3000 = 5.15 (Shen et al. 
2011). After subtracting the pseudo-continuum model 
from the spectrum, we fit a two-Gaussian model to the 
Mg II line and derive a full width at half maximum of 
FWHM = 3670+405 km s~! and a Mg II-based redshift 
of zugrr = 7.611 + 0.004. The detailed spectral fitting 
result is shown in Figure 1. Based on the spectral fitting 
and the Mg II virial estimator of Vestergaard & Osmer 
(2009), we estimate the mass of the central SMBH to 
be Mgg = (1.6 + 0.4) x 10? Mo. The Eddington ratio 
of this SMBH is Lyg/Lgaa = 0.67 + 0.14, which indi- 
cates that the quasar is in a rapidly accreting phase, 
similar to other luminous quasars known at z > 7 (e.g. 
Mortlock et al. 2011; Banados et al. 2018; Wang et al. 
2018; Yang et al. 2020a). Following Wang et al. (2018, 
2020a) and Yang et al. (2020a), all the uncertainties are 
reported as the 16th and 84th percentiles of the distri- 
bution of each quantity, as derived from spectral fitting 
of 100 mock spectra. These mock spectra were created 
by randomly assigning Gaussian noise at each spectral 
pixel, scaled to that pixel’s 1l-o error. All the measure- 
ments and corresponding uncertainties from the spectral 
fitting are listed in Table 1. 

In Figure 1, we also show the individual fittings of the 
100 mock spectra. Note that our best-fit model slightly 
over-estimates the flux in the ranges of 22800-23400 A 
and 24400-25000À. In order to explore the effects of 


the continuum overestimation on Mpy measurements, 
we re-fit the spectrum in K-band only and estimate 
FWHMmen = 4108 + 473 km s^! and Mpy = (1.9 + 
0.3) x 10? Mo, consistent with that derived from the 
global fitting within uncertainties. This K-band spec- 
tral model gives a better fit over this wavelength range, 
but extrapolating it to the J and H bands shows that 
it overestimates the continuum at wavelengths shorter 
than K, suggesting that the quasar could be slightly red- 
dened by dust and cannot be well modeled with a single 
power law. In addition, the fitting could also be affected 
by the difference between the iron template and the iron 
emission from this particular quasar, the possible Mg II 
absorption from outflowing gas and/or intervening metal 
absorbers (see 84), and the large uncertainties in pixels 
at the red edge of the spectrum. Future space-based ob- 
servations of the spectrum beyond 2.5 jum are needed to 
characterize the dust extinction and to better constrain 
the Mpy. In this work, we adopt the best fit from the 
global fitting as our fiducial model. 

'To compare the constraints of seed black hole masses 
and the growth history of the earliest SMBHs from z > 7 
quasars, we plot the growth history of all known z > 7 
quasars with Mpy measurements in Figure 3 by as- 
suming an Eddington accretion with a 1096 radiative 
efficiency and a duty cycle of unity. In this figure, 
all Mpu were measured using the same Mg II virial 
estimator (Vestergaard & Osmer 2009) and thus have 
the same systematic uncertainties (~ 0.5 dex). J0313- 
1806, the highest redshift quasar known and hosting the 
most massive SMBH at z > 7, poses the most challeng- 
ing constraint on the seed black hole mass. Assuming 
Eddington-limited accretion, if the SMBH started its 
growth at redshift z — 15 — 30 (i.e., ~ 400 — 570 Myr 
growth time), it requires a 10* — 10? Mo seed black hole; 
such a seed is inconsistent with being a Population III 
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Table 1. The observational information and 
physical properties of J0313-1806 


R.A. (J2000) 03:13:43.84 
Decl. (J2000) —18:06:36.4 
J 20.92+0.13 
Ks 19.96+0.16 
W1 20.01+0.11 
W2 20.22+0.28 
g, r, z > 26.17, > 25.91, > 24.63% 
Zps1) Ypst > 23.13, > 22.07 
M1450 —26.13+0.05 
Zic] 7.6423+0.0013 
ZMglI 7.611+0.004 
ZCIV 7.523+0.009 
Avmett—[ct] 1072+140 km s^! 
Avciv—Mell 3080+332 km s^! 
FWHMmett 36704405 km s^! 
OX —0.91+0.02 
AL eot (2.7+0.1) x 1078 erg s^! 
Loo (1.4+0.1) x 107" erg s^! 
Msn (1.60.4) x 10? Mo 
Lva/ Lgaa 0.67+0.14 
FWHM cry 312+94 km s^! 
Fici 0.60+0.16 Jy km s^! 
Low (0.80+0.22) x10? Lo 
S228.4GHz 0.45+0.05 mJy 
Lrir (1.5+0.2) x10’? Lo 
Maust ~ 7x107” Mo 
SFR{cn] 40—240 Mo yr! 
SFRrir 22525 Mo yr! 
texp, ALMA 29 min 
texp, FIRE, Longslit 15 min 
texp,Flamingos,JH 30 min 
texp,FIRE,Echelle 362 min 
texp,Flamingos, HK 184 min 
texp,GNIRS 485 min 
texp,NIRES 264 min 


“In these undetected bands, we report the 2c 
limiting magnitudes measured from a 3⁄0 aper- 
ture centered at the quasar position. 


star remnant (e.g. Madau & Rees 2001) or the product of 
dynamical processes in dense star clusters (e.g. Portegies 
Zwart et al. 2004). Instead, direct-collapse black holes 
(DCBH) forming in pre-galactic dark matter halos (e.g. 
Begelman et al. 2006) is the preferred seeding scenario. 
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Q9 E — J0038—1527 at z = 7.03 
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Figure 3. Black hole growth track of z > 7 quasars with 
the assumption of Eddington accretion and a radiative ef- 
ficiency of 0.1 at all times. The curve gives the mass of 
seed black hole required to grow it to the observed mass 
of each SMBH in these quasars. The new quasar, J0313- 
0318, puts the strongest constraint on the seed black hole 
mass, compared to other z > 7 quasars, including J0252- 
0503, J0038-1527, J1243+0100, J1120+0641, J1007+2125, 
and J1342+0928 (see legends, Yang et al. 2019; Wang et al. 
2018; Matsuoka et al. 2019b; Mortlock et al. 2011; Yang et 
al. 2020a; Bañados et al. 2018). The masses of the central 
SMBHs in the previously known quasars are collected from 
Matsuoka et al. (2019b), Wang et al. (2018, 2020a,b), and 
Yang et al. (2020a). The blue, green and yellow shaded re- 
gions define the rough mass ranges of the seed black hole 
produced by Pop III stars, dense star clusters, and direct- 
collapse black holes (DCBH), respectively. 


4. EXTREMELY HIGH VELOCITY OUTFLOW 


Another notable feature of J0313-1806 is that the 
spectrum contains several broad absorption line (BAL) 
features. The BAL features are thought to be produced 
by strong outflows launched from the accretion disk of 
the accreting SMBH. In Figure 4, we show the normal- 
ized spectra in the velocity space of the Si IV, C IV, and 
Mg II emission lines. The composite spectrum (blue line 
in Figure 1) that was used for normalizing the observed 
spectrum was constructed from the spectra of a sample 
of low-redshift quasars with similar relative blueshifts 
between C IV and Mg II lines and similar equivalent 
widths of C IV to those of J0313-1806, following the al- 
gorithm developed by Wang et al. (2020a). It was scaled 
to match the observed spectrum in the observed-frame 
2.0 — 2.2 um. Using the normalized spectrum shown 
in Figure 4, we identified two C IV absorption troughs 
(highlighted orange regions in the middle panel of Figure 
4) at extremely high velocities of (0.171—0.186)c (trough 
A) and (0.109-0.155)c (trough B). We use the “balnic- 
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ity” index (BI; Weymann et al. 1991) to estimate the 
strength of the BALs in the quasar. The measured BIs 
are 1300 km s^! and 5400 km s^! for troughs A and B, 
respectively. Trough B also has a Si IV absorption in 
the corresponding velocity range. The associated Si IV 
absorption for trough A falls into the Gunn-Peterson 
trough and is thus undetectable. Because the potential 
BAL absorptions from Si IV could clobber the proximity 
zone, we can not use this quasar to perform the damping 
wing studies. 

We also considered the alternative explanation that 
the two troughs are Si IV absorption troughs (high- 
lighted orange regions in the top panel of Figure 4) 
with slightly lower velocities (i.e. < 0.1c) We rule 
this explanation out given that there are no associated 
C IV absorption troughs at the corresponding velocities. 
We note that there is a possible weak Mg II absorp- 
tion (highlighted purple regions in the bottom panel 
of Figure 4), which would mean that J0313-1806 is a 
Low-ionization BAL (LoBAL). However, the Mg II ab- 
sorption feature does not satisfy the BI definition (i.e. 
continuously smaller than 0.9 for more than 2000 km 
s71). The absorption could also be affected by a mis- 
match of the iron emission between the composite and 
J0313-1806. Foreground absorption from the interven- 
ing IGM and/or circumgalactic medium (CGM) could 
also contribute to some of the absorption. Future high 
resolution spectroscopic observations are needed to iden- 
tify individual foreground metal absorbers. 

Extremely high velocity (> 0.1 c) outflows are among 
the most promising evidence for active galactic nucleus 
(AGN) feedback, given that their kinetic power might 
be high enough to affect the star formation in quasar 
host galaxies (e.g. Chartas et al. 2009). Such outflows 
are also very rare phenomena at lower redshifts; for ex- 
ample, Rodríguez Hidalgo et al. (2020) identified only 
40 (~ 0.6%) extremely high velocity outflow systems 
from a parent sample of — 6700 luminous quasars at 
z ~ 2— 5. Thus, the discovery of relativistic outflows 
in J0313-1806 and J0038-1527 (z = 7.03; Wang et al. 
2018) among a sample of only 8 z > 7 quasars indicates 
that relativistic outflows are more common at the high- 
est redshifts, suggesting an evolution with redshift of 
the quasar outflow phenomenon. In addition, the C IV 
broad emission line has a substantial blueshift, another 
signature of radiation driven outflows, with a velocity 
of 3080 km s^! relative to the Mg II line and 4152 km 
s^! relative to the [C II] line. The substantial blueshift 
of the C IV broad emission line, consistent with what is 
expected from the relation between the outflow velocity 
and the C IV blueshift (Rankine et al. 2020), is among 
the most extreme cases observed at lower redshifts. Such 
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Figure 4. The normalized spectrum of J0313-1806. The 
solid and dashed lines represent 10096 and 9096 of the nor- 
malized spectrum, respectively. The top panel shows the 
outflow velocity of the Si IV line, the middle panel denotes 
the outflow velocity of the C IV line, and the bottom panel 
represents the outflow velocity of the Mg II line. We interpret 
the two most obvious troughs (orange regions) as extremely 
high velocity C IV outflow systems with trough A having v — 
(0.171 — 0.186)c and trough B having v — (0.109 — 0.155)c. 
There is a potential weak Mg II absorption trough (the high- 
lighted purple region in the bottom panel), but it does not 
satisfy the BI definition. 


large blueshifts of C IV emission have also commonly 
been observed in other luminous z > 6.5 quasars (e.g. 
Meyer et al. 2019; Mazzucchelli et al. 2017; Wang et 
al. 2020b), but have not been commonly found in lower 
redshift quasars. These observations indicate that the 
earliest SMBHs are in a fast-growing phase (i.e. with 
a high Eddington ratio) and that quasar outflows could 
play a crucial role in regulating the growth of the earliest 
SMBHs and their hosts. 


5. QUASAR HOST GALAXY 


The dust continuum around the redshifted [C II] emis- 
sion from the quasar host galaxy is significantly de- 
tected by ALMA (Figure 2). We measure the contin- 
uum flux S228.4 GHz = 0.45 + 0.05 mJy from a two- 
dimensional Gaussian fit to the collapsed continuum 
map using the CASA task imfit. The continuum emis- 
sion is marginally resolved with a de-convolved size of 
0"62 x 0729 or a physical size of 3.12 x 1.44 kpc, similar 
to the sizes of other quasar host galaxies at z = 6 (e.g. 
Venemans et al. 2018). Under the standard assumption 
of an optically thin greybody dust spectral energy dis- 
tribution with a dust temperature of Taust = 47 K and 
an emissivity index of 8 = 1.6 (Beelen et al. 2006), af- 
ter considering the effect of the cosmic microwave back- 
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ground on the dust emission, we estimate a total infrared 
luminosity of Lr = (1.50.2) x 10!? Lo. We estimate 
the star formation rate (SFR) of the quasar host galaxy 
to be SFR=225+25 Mo yr-! based on the scaling rela- 
tion between the SFR and Lrır (Murphy et al. 2011). 
With the same set of assumptions, we estimate a total 
dust mass of ~ 7 x 10’ Mọ following Venemans et al. 
(2018). 

The [C II] line is also detected (with peak pixel at 
> 4c) in our ALMA observations as shown in Figure 
2. A single Gaussian fit to the [C IT] line gives zor = 
7.6423 + 0.0013 and FWHM cry = 312+ 94 km sl. In- 
stead of measuring the total [C II] line flux from the one- 
D spectrum, we measure it from the two-dimensional 
integrated line map using imfit following Decarli et 
al. (2018), and derive Fic = 0.60 + 0.16 Jy km s^ !. 
This corresponds to a [C II] luminosity of Lior = 
(0.80+0.22) x 10? Lo and a SFRjorj = 40—240 Mo yr^! 
by adopting the empirical relation from De Looze et al. 
(2014), consistent with the SFR derived from the total 
infrared luminosity. The properties of the quasar host 
galaxy are comparable to that of the other two quasar 
host galaxies known at z > 7.5 (Venemans et al. 2017; 
Yang et al. 2020a). Future high-resolution observations 
of the star formation distribution with ALMA and ob- 
servations of the kinematics of extended [O III] emis- 
sion with the Near-Infrared Spectrograph (NIRSpec) on 
James Webb Space Telescope (JWST) will be necessary 
to carry out more detailed investigations of the dynami- 
cal mass of the host galaxy and the impact of the quasar 
outflow on star formation in the host galaxy. 


6. SUMMARY 


In this letter, we report the discovery of a luminous 
quasar, J0313—-1806, at redshift z = 7.642. It is the most 
distant quasar yet identified. J0313-1806 has a bolomet- 
ric luminosity of Lyo; = (1.4 + 0.1) x 104” erg s^! and 
hosts a SMBH with a mass of (1.6 + 0.4) x 10? Mo, 
accreting at an Eddington ratio of 0.67 + 0.14. The 
existence of such a SMBH just ~ 670 million years af- 
ter the Big Bang puts strong constraints on the forma- 
tion models of seed black holes. The quasar’s rest-frame 
UV spectrum exhibits broad absorption troughs from 
extremely high-velocity outflows. These outflows have 
a maximum velocity up to ~ 20% of the speed of light. 
We also detect strong dust emission and [C II] line emis- 
sion from the host galaxy in ALMA data. The ALMA 
observations suggest that J0313-1806 is hosted by an 
intensely star-forming galaxy with a star formation rate 
of ~ 200 Me yr-!. The continuum observations in- 
dicate that substantial dust (~ 7 x 107 Mo) was al- 
ready built up in the quasar host galaxy. The relativistic 


quasar outflow and the fast SMBH growth phase, com- 
bined with the intense star forming activity in the host 
galaxy, suggest that J0313-1806 is an ideal target for 
investigating the assembly of the earliest SMBHs and 
their massive host galaxies with future high resolution 
ALMA and JWST NIRSpec/IFU observations. 
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